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Abstract 

We study the frictional sliding of two ideally incommensurate surfaces with a third 
incommensurate sheet - a sort of extended lubricant - in between. When the mu- 
tual ratios of the three periodicities in this sandwich geometry are chosen to be the 
golden mean (/> = (! + \/5)/2, this system is believed to be statically pinned for 
any choice of system parameters. In the present study we overcome this pinning 
and force the two "substrates" to slide with a mutual velocity Vext, analyzing the 
resulting frictional dynamics. An unexpected feature is an asymmetry of the rela- 
tive sliding velocity of the intermediate lubricating sheet relative to the two sub- 
strates. Strikingly, the velocity asymmetry takes an exactly quantized value which 
is uniquely determined by the incommensurability ratio, and absolutely insensitive 
to all other parameters. The reason for quantization of the velocity asymmetry will 
be addressed. This behavior is compared and contrasted to the corresponding one 
obtained for a representative cubic irrational, the spiral mean lo. 



1 Introduction 



The intriguing concept of incommensurability has recently found relevant and 
practical applications in the discipline of nanotribology, i.e., the science of 
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friction, lubrication and wear at the nanoscale, in particular for the possibility 
of giving rise to superlubricity [1]. When two hard crystalline surfaces with 
incommensurate lattice spacings (or commensurate but not perfectly aligned) 
are brought into contact, then the minimal force required to achieve sliding, 
known as the static frictional force Fs, vanishes, provided the two substrates 
are stiff enough. This remarkable conclusion can be drawn, for example, in 
the context of the Frenkel-Kontorova (FK) model [2,3]. In real situations how- 
ever such a case of "dry" friction is exceptional. A physical contact between 
two solids is generally mediated by so-called "third bodies", which act like 
a lubricant film. The sliding interface is thus characterized by three distinct 
length scales: the periods of the bottom and top substrates, and the period 
of the lubricant structure. Upon relative sliding of the two substrates, physi- 
cal intuition would suggest that this confined lubricant should move with the 
average of the two substrate velocities. At the nanoscopic level, however, in- 
commensuration between the atomic corrugations of the solid substrates could 
change this expectation, with possible tribologically important implications of 
an irregular distribution of the lubricant velocity on the wear of the sliding 
surfaces. 



In order to study the role of incommensurability among the three interface 
inherent lengths on the sliding dynamics, we consider, as in Ref. [4], a sim- 
plified one-dimensional generalized FK model consisting of two moving rigid 
sinusoidal substrates, of spatial periodicity a_|_ and a_, and a chain of harmon- 
ically interacting particles, of equilibrium length a^, mimicking the lubricant 
and sandwiched between the two substrates, as shown in Fig. 1. 



In the standard FK model (single substrate potential, of period a_|_, plus har- 
monic chain), the static behavior depends strongly on the value of r = a+/ao: 
a rational r produces pinning (a finite force is necessary to induce lubricant 
translation), while an irrational r is associated to a phase transition (the 
Aubry transition) between a pinned phase, for small spring stiffness K, to an 
unpinned regime, for large iiT [5]. A clear critical behavior was demonstrated 
only for quadratic algebraic irrational numbers, among which the Aubry crit- 
ical Kc for depinning is smallest for r = = (1 -|- •\/5)/2, the golden mean 
(GM) [6]. In this sense, therefore, the GM represents the "most irrational" 
incommensuration, with the largest unpinned phase. 



We will show that, on the contrary, the dynamics of our generalized FK model 
shows a surprising result, the Golden Mean case leading to a strong dynamical 
commensur ation . 
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2 Model 



Consider the system sketched in Fig. 1. The equation of motion of the i-th 
lubricant particle is: 

mxi — — - [F+ sin — v+t) + F_ sin k-{xi — V-t)] 

+ K{Xi+i + Xi^i - 2Xi) - 2j{Xi - Va) , (1) 

where m is its mass, K is the chain spring constant, k± — (27r)/a± are the 
wave- vector periodicities associated to the two sinusoidal substrates of ampli- 
tude U±, moving at velocities v±. Here Va = {v^+V-)/2 is the average speed of 
the substrates, and 7 is a viscous friction constant accounting phenomenolog- 
ically for various sources of dissipation in the substrates (electrons, phonons, 
etc.). (The last term in Eq. 1 originates from two frictional contributions of 
the form —^{xi — v^) — j{xi — f-)-) F± = k±U± are the amplitudes of the 
forces due to the corrugation of the substrates, and we set F_/ Fj^ = 1 as the 
least biased choice. The spring equilibrium length ao (not entering explicitly 
the equation of motion (1)) enters through the boundary conditions, which 
we take to be periodic (PBC), xnq+i — xi -\- NqUo- The traditional imple- 
mentation of the PBC [4,7] is realized by approximating a_/a+ and a+/ao 
(via, for instance, a continued fraction expansion [8]) with suitable rational 
numbers N_/N^ and N^/No numerically close to r (e.g. successive integers 
N_, N_^, Nq in the Fibonacci sequence approximate r = (f)). Convergence of all 
physical quantities is then checked in the limit of large N± and A^o, where full 
incommensurability is restored. In the present case, however, the boundary 
condition enters only in the K-term, and there is no real need to approximate 
r: we can use the machine-precision value of r. This alternative procedure in- 
troduces, at the boundary, a small phase "twist" to the purely local substrate 
forces, but this boundary error can be made arbitrarily small by choosing a 
suitably large number Nq of particles. All the results presented below coincide 
for either implementation of the PBC. 

The general behavior of the driven system in Eq. (1) depends crucially on the 
nature of the incommensurability of the substrates and the chain. In order 
to make a comparison with previous tribological studies [4,7], as well as to 
previous studies on generalized FK models [9], we restrict our present anal- 
ysis to cases where the ratio a_/a+ of the substrate periodicities and the 
ratio a+/ao, fixing the equilibrium length of the lubricant chain, are equal: 
r = a_/a+ = a+/ao > 1- We take a+ = 1, m = 1, and F_|. = 1 as basic 
dimensionless units. 

The equations of motion are integrated using a standard fourth-order Runge- 
Kutta algorithm. The system is initialized with the chain particles placed at 
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rest at uniform separation oq. After relaxing the starting configuration and 
selecting a reference frame in which the substrate is at rest (f+ = 0), 
the top substrate U_ begins to move at the constant velocity v_ = V^xt- 
After an initial transient, the system reaches a "dynamical stationary state" , 
characterized by generally fluctuating but significantly (on average) constant 
drift velocity of all particles. 

Deferring a more general discussion to later work, we address here the con- 
trasting behavior of two different incommensurabilities r, namely, the golden 
mean (GM) and the spiral mean (SM) , belonging to distinct classes of 
irrational numbers (quadratic and cubic, respectively). The very different be- 
havior of golden and spiral will demonstrate how strongly the precise value of 
incommensurability dramatically influences the tribological behavior of inter- 
face sliding. 

Statically [9,7], it was found that when the length scale ratio r equals the SM 
there exists a value of the harmonic interaction strength K, above which the 
static friction force Fg is zero. On the contrary, when r is equal to the GM, 
then Fg is nonzero for any value of K: the system is always pinned. 

Dynamically, the same model was recently considered with the two substrates 
forced to slide relative to one another by an external force applied through an 
elastic spring [4] . It was found that the GM and the SM behaved very differ- 
ently, showing that there is no quantitative or even qualitative uniformity of 
behavior in incommensurate kinetic friction, and that certain incommensura- 
bilities favor sliding systematically better than others. 

In the present study of dynamic sliding we focus on the motion of the lubri- 
cant film and analyze the corresponding particle motion. The sliding of one 
substrate relative to the other is rigidly forced, in the form of a fixed relative 
velocity (no elastic pulling, as in Ref. [4]). 



3 Results and discussion 

Figure 2 shows the time-averaged lubricant chain center of mass (CM) velocity 
Vcm: after stationarity has been reached, as a function of the chain stiffness 
K for the GM and SM values of r. The value of 7 is here fixed to be 7 = 0.1 

^ The golden mean cj) = {^/5 + l)/2 ~ 1.6180 is the irrational number satisfying the 
algebraic quadratic equation 0^—0—1 = 0; its rational approximants [8] are provided 
by the Fibonacci sequence, Fn-^-i = Fn + Fn-i, with Fq = Fi = 1: Fn+i/Fn 4>. 
The spiral mean (u ~ 1.3247) belongs to the class of cubic irrationals, satisfying the 
equation — w — 1 = 0; its rational approximants can be generated by the recursion 
relation Gn+i = Gn-i + Gn-2 with G_2 = Go = 1, G-i = 0: Gn+i/Gn —>■ uj 
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(under-damped regime), and the external velocity is Vext = 0.1. The first 
obvious observation is that only in the limit of large K does the lubricant 
move at the mean speed Va = Vext/'^'- generally, the CM velocity Vcm is a 
complicated function oi K. It encompasses regimes of continuous evolution, 
but also remarkably and unexpectedly flat plateaus. 

The most striking feature is indeed the presence of perfectly flat Vcm/yext 
plateaus, whose precise value is remarkably independent not only of but 
also of 7, Vext-, and even (we checked) of F^/F^. Super flcially, the two in- 
commensurability considered, GM and SM, appear to share the same general 
features. Going deeper into the details of the particle motion, the two cases ap- 
pear to be somewhat different. In order to clarify these differences, we plot in 
Figs. 3 and 4 the typical time-evolution of the velocity of a particle and of Vcm-i 
for a value of K inside a plateau: K — 1. For the GM case (r = 0), the motion 
of each single particle is perfectly time-periodic with period Vext T c:^ 2.62 (see 
Fig. 3, in particular the Fourier spectrum panel). Like Km/Kxt, also VextT 
is found to be independent of all parameters. The motion of every particle in 
the chain is simply characterized by a phase- variable 6i, so that we can write, 
in the stationary state, Vi{t) = g{ujt — 6i) where g{z) is a periodic function in 
[0, 27r], and u — 2t^/T. On the contrary, the motion is not periodic in the SM 
clear from the Fourier spectrum panel in Fig. 4. 

The key to the exact quantization, periodicity of motion, and universality of 
Vcm/Vext and Vext T for the GM case, is what we might call a form of dynamical 
commensurate sliding, which, surprisingly, the GM incommensurability shares 
with all the rational values of r. We now illustrate this effect. Consider travel- 
ing with a given center of mass velocity < Vcm < Vext- Correspondingly, let 
/+ = T-^ = Vcm/a+ and /_ = Tz^ = {Vext - Vcm) / d- be the average frequency 
of encounter with the substrate corrugation minima. It turns out that, for the 
GM case, the Vcm is perfectly tuned so as to make the two frequencies /+ and 
/_ mutually commensurate, indeed perfectly equal: 

/+ = /_ (Golden Mean matching) , (2) 



so that the particles can "satisfy" both substrates by a concerted motion with 
period T = f^^ = a^/Vcm- It is a matter of simple algebra to realize that 
Eq. 2 is fulflUed when the average lubricant speed Vcm is given by: 

Vcm _ Q+ _ 1 /o\ 

Vext a+ + a_ l + ' 



Correspondingly, the period of the motion is given by VextT = a+(l + 0), and 
the fundamental frequency of the Fourier spectrum by /+ = Vextd^^/il + 4>)- 
Eq. 3 describes precisely the value of the plateau in Vcm,/Vext show in Fig. 2 for 
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the GM case. Similarly, the value of /+ is precisely the fundamental frequency 
shown in the Fourier panel of Fig. 3. The GM plateau in Km/Ki^t is remarkably 
robust and wide (indeed, its width as a function of K increases for decreasing 
Vext)'- we have found no other irrational with a similarly wide plateau, which 
definitely competes with those obtained for rational values of r (not shown) . In 
the present three-slider dynamical context, therefore, the GM ratio represents 
the most ^^dynamically commensurate irrational". That is just the opposite of 
the statics of the standard FK model [2,3] where GM incommensurability is 
the hardest one to pin. 

The SM case is quite different. The plateau is not just smaller: it is not quan- 
tized by any frequency matching condition as in Eq. (2). Correspondingly, the 
particle motion in the chain is definitely not periodic, as confirmed by the 
Fourier spectrum panel in Fig. 4. 



4 Conclusions 

When two periodic substrates are forced to slide relative to one another but 
there is an interposed "lubricant" layer in between, the velocity of the inter- 
mediate layer may display an interesting behavior. In general, this velocity 
will not be the average between that of the two substrates. Moreover, in some 
cases at least, the lubricant layer velocity may be quantized to an exact value, 
depending only on relative commensur abilities. That occurs not only if the 
two substrates and the layer are mutually commensurate (rational values of 
r = a_/a+ = a+/ao), but strikingly also when they have certain types of in- 
commensurabilities, for instance their periods are related through the golden 
mean, r = 0. On the other hand, that certainly does not apply to all types 
of incommensurability, and we found no simple matching quantization and a 
radically different behavior (aperiodic motion) in the case of the spiral mean, 
r — uj. Work is under way to fully explore this interesting scenario. 



Acknowledgments 

We are grateful to O.M Braun for collaboration, stimulation, and illumi- 
nating discussions. This work was partly supported by MIUR COFIN No. 
2003028141, MIUR COFIN No. 2004023199, by FIRB RBAU017S8R oper- 
ated by INFM-CNR, by MIUR FIRB RBAU017S8 R004, and by INFM-CNR 
(Iniziativa trasversale calcolo parallelo). A.V. was supported in part by PRRI- 
ITT (Regione Emilia Romagna), Net-Lab "Surface & Coatings for Advanced 
Mechanics and Nanomechanics" (SUP&RMAN). 



6 



References 



[1] M. Dienwiebel, G.S. Verhoeven, N. Pradeep, J.W.M. Prenken, J. A. Heimberg, 
and H.W. Zandbergen, Phys. Rev. Lett. 92, 126101 (2004). 

[2] S. Aubry, Physica D 7, 240 (1983). 

[3] S. Aubry and P.Y. Lc Dacron, Physica D 8, 381 (1983). 

[4] O.M. Braun, A. Vanossi, and E. Tosatti, Phys. Rev. Lett. 95, 026102 (2005) 

[5] M. Peyrard and S. Aubry, J. Phys. C: Solid State Phys. 16, 1593 (1983). 

[6] J. Greene, J. Math. Phys. 20, 1183 (1978). 

[7] A. Vanossi, J. Roder, A.R. Bishop, and V. Bortolani, Phys. Rev. E 63, 017203 
(2001). 

[8] A.Ya. Khinchin, Continued Fractions (Dover, New York, 1997). 

[9] T.S. van Erp, A. FasoHno, O. Radulescu, and T. Janssen, Phys. Rev. B 60, 
6522 (1999). 



7 





Fig. 1. Schematic drawing of the model with three characteristic length scales. Here 
oo is to be interpreted as the equilibrium distance of the harmonic chain (lubricant 
film). a± are the spatial periods of the two sinusoidal substrate potentials, assumed 
to be rigid. 
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Fig. 2. Normalized velocity of the center of function of the 

chain stiffness K, for golden mean (GM) and spiral mean (SM) incommensurability. 
Here 7 = 0.1 and Vext = 0.1. Note the logarithmic scale in the abscissa. Observe 
that Vcm/yext approaches the expected value 1/2 only for large values of K (hard 
lubricant). The value of the wide GM plateau is exactly tuned by a dynamical 
matching condition which results in Vcm/yext = 1/(1 + 0) ~ 0.381966 (see Eq. 3). 
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Fig. 3. Typical time-evolution of the velocity of a particle Vi, and of the cen- 
ter of mass Verm for a value of K inside the GM plateau {K = 1). T indi- 
cates the time periodicity of the sliding dynamics at stationarity, exactly given 
by T = a+(l + <j))/Vext ~ 26.180 (see Sec. 3). The lower panel displays the extraor- 
dinarily sharp peaks of the corresponding Fourier transform of Vi, showing a single 
frequency, exactly located (see Sec. 3) at /+ = T^^ = VextO'+ /(I + ^ 0.0381966 
(plus higher harmonics). 
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Fig. 4. Same as Fig. 3, but for the SM incommensurability, still for a value of K 
inside a plateau {K = 1). At variance with Fig. 3, no periodic dynamics is here 
observed, as the noisy corresponding Fourier spectrum (lower panel) gives further 
evidence of. 



11 



